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Abstract 

We generalize effective energy variational techniques to study appropriately 
Qh! quantized solitonic field configurations. Our approach rests on collective quan- 

r~| , tization ideas and is specifically designed for the numerical evaluation of soliton 

parameters. We employ this method to obtain the one- loop quantum corrections 
K> , to the soliton mass and form factor. Special attention is given to the regularization 

H , of the physical observables in the solitonic sector of the theory. The numerical im- 

plementation of the method is demonstrated for a simple one-dimensional scalar 
field example. 
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1 Introduction 

The dilemma between our classical visualization of a soliton as a compact object and the 
plane-wave nature of its energy and momentum quantum eigenstates is as common as 
the quantum description of any composite particle. Not surprisingly, this problem has 
been approached with diversity of techniques'^"^' such as functional integral, canonical 
quantization, semi-classical and variational procedures. The methods that have already 
been developed are able to provide reasonable qualitative understanding of the physics 
behind a "quantized soliton" but, as a rule, they are analytically tractable only for 
the simplest models, mostly in one space dimension. Whenever it comes to definite 
predictions for physical quantities, many ambiguities appear'^'^': consistent ultra-violet 
cutoff methods, vacuum energy subtraction, and boundary conditions at infinity being 
only some of the examples. These questions must be settled if we want the method to 
yield a definite finite result. 

The present paper is intended to complement the earlier works on soliton quantization 
in the following three aspects. First, the approach that we develop naturally unifies the 
classical image of a soliton as a solitary wave and the fundamental quantum-mechanical 
principles. It makes a simple connection with such more formal views on this problem as 
the collective quantization methods'^'^' and the form-factor, "Kerman-Klein" , approach 
applied to solitons in the Goldstone and Jackiw paper'^l This is demonstrated in Sec- 
tions 2-4, where we review the basics of collective quantization (Sec. 2), then introduce 
our variational scheme (Sec. 3), and draw parallels to the Kerman-Klein method for 
calculating the soliton form- factor (Sec. 4). 

Second, during this work we have had in mind to construct a calculational method 
for soliton parameters which can be practically implemented as a computer algorithm. 
It then could be applied to analyze the quantum effects in theories that are too compli- 
cated for analytical solution, but from qualitative considerations'^' may possibly possess 
non-trivial stable solitonic states. Much progress in this direction has already been 
achieved'^'^'^l in the last couple of years. However, the developed numerical methods for 
calculation of quantum corrections would not be fully satisfactory without understand- 
ing of the special role played in the system description by the cyclic variables, such as the 
soliton position in space. We discuss the relative significance of the modifications needed 
to incorporate the cyclic variables (Sees. 3,6.2) and show on a simple example of a one 
dimensional kink (Sees. 6.1,6.2) how to calculate the one-loop quantum corrections to 
the physical parameters characterizing the soliton in its true, delocalized, ground state. 

Third, we discuss theory regularization in the solitonic sector remembering that 
renormalization conditions are conventionally imposed in the perturbative sector of the 



theory (Sec. 5). In dealing with this subtle issue, we prefer to apply physically realizable 
regularization methods rather than formal manipulations with divergent expressions. 
Our final formulas involve only finite quantities and convergent integrals, suitable for 
the numerical computations of Section 6. A summary and conclusions are given at the 
end of the paper (Sec. 7). 

2 Collective quantization 

We assume that the theory has an absolutely stable solitonic state. That is, the decay 
of this soliton into plane-wave excitations over the usual vacuum is forbidden by con- 
servation of some charge Q, which could be a topological charge or a fermion number. 
Different values of the charge Q split the overall Hilbert space of the theory into sepa- 
rate non-mixing sectors with their own ground states, which would be the true vacuum 
for the "perturbative" sector and the soliton at rest for the "solitonic" sector. Due to 
translational invariance of the theory, the total momentum operator P commutes with 
the Hamiltonian, 

and the lowest energy state in the solitonic sector should also be an eigenstate of the 
momentum operator, corresponding to the total momentum zero. By the Heisenberg 
uncertainty principle, the position of the soliton in space is completely undefined. The 
techniques to handle this situation are known'^^^1 . 

To be specific, we illustrate the formalism on the simplest example of a real scalar 
field (f) in one space dimension with the Lagrangian density 

This theory may support a stable topological soliton if the global minimum of the po- 
tential U is not unique: 

3 cTi, a2 : U{ai) = U{a2) = min W , ai ^ a2 

(we assume that the minima cxi and (J2 are connected by a discrete symmetry of the 



Lagrangian (|2.1| ), e.g. 0^—0.) Let (j{x) be an arbitrary function such that 



lim a(x) = cxi , lim a(x) = 02 ■ (2-2) 



Let also {^A:(a^)} be a set of orthonormal functions, 

dx rik{x)rii{x) = 6ki 



that are orthogonal to a"'(x), 

dx cr'(x) rik{x) = , 

and, together with cr'{x), form a complete set. 

Following the Christ and Lee method^, we describe the solitonic sector of the theory 
trading the field 0(a;) for an equivalent set of commuting dynamical variables {X, gi, q2, ...} 
as 



{X 



g V(a; - X) + ^ qkVk{x - X) . 



(2.3) 



A;=l 



Given a field configuration 0(x), the corresponding variables {X, gi,g2, •••} may be de- 
termined from the equations 



dx a'{x — X)(j){x) = , 
qk= dx rik{x - X) U(x) - 5'"V(x - X) 



(2.4) 
(2.5) 



Let {P, vTi, 71-2, ...} be the operators of the canonical momenta conjugate to {X, qi,q2, ...}, 
so that 



[P,X] 



Ffc, qi\ 



-i5^ 



kl 



with all other possible commutators being zero. The original Hamiltonian, 

1 



with 



H[t^{x),<P{x)] = - I dx {iT{x)y + V[(f)] 



v[<p] 



dx 



-{<P'{x)f + g-'U{g<P{x)) 



(2.6) 



(2.7) 



when being re-expressed in terms of the new variables and their conjugate momenta, 
reads'^'-^^l : 



-| oo 

H{P,7r,q) = Y.^l + Viq)+9'h{P,n,q) 
fc=i 



(2.^ 



In this expression h{P, vr, q) is an 0{g'^) correction to the kinetic term that is quadratic 
in the canonical momenta P or vr, but has somewhat intricate g-dependence''*', and V{q) 
is simply the potential term ( |2.7|) written in terms of the new variables. 

The Hamiltonian ( p.8| ) does not depend on the "collective coordinate" X to any order 
in the coupling g reflecting the system invariance to translation of the field 0(x) in space 



as a whole. Its conjugate momentum P is, therefore, a conserved quantity: 

[P,H] = 0. 

It is possible to showt'^'^^' that the operator P indeed represents the total momentum 
of the systemQ. The essence of the collective quantization is to describe the original 
theory restricted to the subspace of the quantum states \Qp) that carry a definite total 
momentum P, for example 

p\no) = . 

The zero value of the momentum can always be achieved by a proper choice of the 
Lorentz frameQ. In the P = sector the full Hamiltonian ( |2.8| ) reduces to 

H'{fr,q) = -Y.^l + V{q) + 0{g'), (2.9) 

k=l 

and the remaining dynamical variables may be handled by conventional perturbative or 
semi-classical methods. 

We would like to conclude this section by the following remark: In the original Christ 
and Lee method the profile function g^^a{x) in eq. (|2.3|) is required to satisfy the classical 
field equations of motion. In the present discussion, a{x) is an arbitrary function with 
the proper asymptotic behavior, eq. ( p.2|) , yielding a non-degenerate transformation 
(t){x) -^ {X, gi, g2, •••} • We exploit this freedom in the choice of a soon. 

3 Reduced effective energy 

For the solitonic sector of the theory and a given function 0"(x) in the canonical trans- 
formation fBTSD, we define 



{j[dq\ e*'^ wl) J 
r'[gi(t),g2(t),...] = min { -W' [J] + Y.fdtJ,{t)qu{t)\ (3.2) 



^ Let us provide a simple illustration of this statement observing that the commutation relations 
[P,qk]=0 and [Pj{X)] = -if'{X) V function /(x) 
being applied to eq. ( |2.3| ) produce immediately the formula for the total momentum commutator: 

^ Quantization in moving frames may also be consideredl^^^ . 



where S'[q] = S'[qi{t),q2{t), ...] in eq. (|3.1|) is the classical action of our system in the 
P = sector corresponding to the "reduced" Hamiltonian (|2.9| ), and the functional mea- 
sure [dq] = rifcli [d-Qk] absorbs all the factors that may have appeared in the functional 



integral when the canonical momenta tt^ were integrated out. Eqs. ( p.l|) and ( |3.2| ) are 
reminiscent of the standard definitions of the generating functional iy[J(x)] and the 
effective action r[0o(a^)], when the source terms would have the form J d'^x J{x)(f){x). 
Here, we couple the sources Jk to the complicated non-/mear functionals gjt[0(x)]. We 
refer the quantity V as the "reduced effective action" . 

Let us consider the special case when the arguments of V are constant functions 
qi,q2,---, and let the time integrals in eqs. (|3.1|- p.2| ) be taken over a large but finite 
interval < t < T. Then 

r'[gi(t), q2{t), ...] = E'^ffiqi, g2, •••) T , 

where the reduced effective energy Ecs allows a nice physical interpretationt^^l . Namely, 
it is the minimum of the Hamiltonian expectation value among the quantum states \Q') 
for which the quantum expectation values of the operators {qk} are given by E'^f^ argu- 
ments, {qk}: 

{n'\H'\n') 

i^efrfa,g2,-)=„,min^^^ . (3.3) 

We primed the ket- vectors 1^2') in order to emphasize that these states belong to the 
Hilbert space associated with the reduced Hamiltonian H'. It is different from the space 
of the states, {|f2)}, of the original theory described by H because the latter carry an 
additional quantum number - the total momentum P. However, by construction 

H'in,q) = HiP,n,q)\^^^ , 

and one can rewrite eq. ( |3.3|) in terms of the objects of the original theory: 

F' (n n )- miu ^^1^'^) 

(n|qfc|n)=9fe,p|n>=o 

Provided that we are able to compute the reduced effective energy E'^f^, the mass of the 
soliton in its ground state, M, may be determined as 

M= mm E'^s{qi,q2,...)= min \' '/ . 3.4 

91,92, ■■• ^°^^ ^ ' in) such that (QIQ) 

p\n)=o 

Now we proceed to the actual computation of E'^f^. Our first crucial step is to trade 
the discrete infinite set of variables {gi, q2, ■■■} for a more physically intuitive continuous 

6 



"shape" function. As eq. (p^ ) suggests, the mass of the sohton, in principle, could be 
calculated by searching for the minimum of E'^f^ over the range of parameters {gi, q2, ■■■} 
that are defined by means of the decomposition ( p.3| ) with an arbitrarily chosen but fixed 
function cr(a;). However, we find it more straightforward to always calculate E'^^ at zero 
expectation values of all qk and vary the function a{x) itself looking for the soliton mass 
as 



M = rmnE'^s[a{x)] 

a{x) 



(3.5) 



In this equation, E'^Q[a{x)] is i?^fj(0,0, ...), by our previous notations, where the coordi- 
nates {qk} are defined using the given "shape" function a{x). As shown in Section 4, the 
configuration a{x) that minimizes the right hand side of eq. ( p.5|) is closely connected 
to the soliton form factor. 



The semi-classical expansion for the reduced effective energy 
E'^^[a{x)] = E,,[a{x)] + E'cJa{x)] + 0{g) 



(3.6) 



starts from an 0{g ^) term - the classical energy associated with the field configuration 



[x)=g 'a[x) 



E, 



y[a] = V[g-^a] = -( dx 

U J 



^-ia'ix)f+Uiaix)) 



(3.7) 



),'?] 



The leading quantum correction, j; In det f — ss 

to 0{g) terms neglected in the one- loop approximation) as the "Casimir energy" 



g=0 



0((7°), can be written l-^''! (up 



\k=l / regularized 



(3.8) 



where {uJk}k=i,2,... by definition are the frequencies of small oscillations in the reduced 
system "stabilized" at g = by an "external source". That is, {uj'j^} are the normal 
frequencies of the Hamiltonian 



H'in,q) = H'in,q)-j:\^ 

fc=i \ c'gfc 



Qk 



9=0; 



1 OO OO / QT/ 



qk 



(3.9) 
(3.10) 



g=Oy 



We claim that, with the accuracy including at least 0{g), the frequencies {uj'^.} are 
given by (the square roots of) the non-zero eigen-values in the following Sturm-Liouville 



problem: 



dx"^ 



+ U"{a) 



rin{x) + J dy 
[Jdzf^{z)a'{z)] 



j'A^yiy) + (y'{x)]'M 



jdz{a'{z)f 



Jdzia'iz)) 
2 cr'ix)a'{y) \ 7]n{y) 



(3.11) 



where 



Ja[x) = -g 



SH 



50(x) 



cr"{x) — U' (cr(x)) 



(3.12) 



=9 o- 



One of the several methods to prove this statement consists in modifying the kinetic 
term in the "stabilized" reduced Hamiltonian (|3.1CI|) by addition of the P-dependent 
piece g^h{P,TT,q) that appears in eq. ( p. 81 ). Then we transform the canonical variables 
back to the continuous fields (p{x) and n{x), when eq. ( p.5|) comes helpful. As the result 
of this procedure, 



H'{7T,q) 



J27rl + V{q)+g'h{P,n,q) 



fc=i 




Qk 



(3.13) 



q=Oj 



H[7i{x) 



{x)] + - I dxj^{x)(f){x + X[(t)]) = H[7T{x),(f){x)] , (3.14) 



where H[tt{x), (f){x)] is the Hamiltonian of the original theory ( p.6| ) and the second term 
in eq. ( |3.14]) is the transform of the stabilizing term from the previous expression. The 
non-linear, non-local functional X[0(x)] in eq. ( |3.14| ) is explicitly defined by eq. ( |2.4| ) and 
the source jo-(x) is precisely the one that appears in eq. (|3.12|) . It is not hard to argue 
that the modification ( |3.13D does not affect the spectrum of small oscillations about the 
classical equilibrium point of this system more than 0((7^), although it implies an addi- 
tional degree of freedom associated with the momentum P and the oscillation spectrum 
of H[k{x),(J){x)] should acquire an additional eigen-mode. However, the new canonical 
variable , "X" , does not explicitly enter H as it was not present in the previous expres- 
sion (p.l3| ), leaving the Hamiltonian (|3.14| ) still being invariant under 0(x) ^ (j){x + a) , 
and the eigen-frequency of the new, translational, mode being identically zero for any 
configuration a{x). Finally, the oscillation spectrum of H is determined by the eigen- 
value problem 



(^n)Sn(a;) 



dy 



S^H 



(50(x)(50(y) 
that after some algebra leads to eq. (|3.11| ) above. 



Vn{y) 



=g o"y 



To summarize, the one-loop contribution to the effective energy in the P = sector 
is given by the sum0 of ^^^^[cr] that are determined by the Schroedinger-type equa- 
tion ( p.ll| ) with a non-local but separable potential. This sum is terribly divergent and 
should be regularized by subtracting the Casimir energy of the trivial vacuum and proper 
counterterms specified by renormalization conditions in the perturbative sector. We per- 
form this regularization in Section 5. In the absence of the non-local terms in eq. ( |3.11| ), 



the sum of the corresponding eigen- modes, including ujq, would yield the one- loop con- 
tribution to the conventional effective energy associated with the spatially "nailed down" 
field configuration 0(x) = g~^(j{x) . The non-local terms vanish for a{x) = crci(x) = g 0ci(a^) 
that solves the classical Euler-Lagrange equation, 

^-W'(a,i) = 0, (3.15) 

because then j^r^iix) =0 by its definition (|3.12| ). These terms guarantee that eq. (|3.11|) 



allows a normalizable zero-frequency solution t]o{x) occr'(x) even when cr(x) differs from 
the classical configuration. However, the effect of the non-local terms on the frequencies 
with n 7^ is only quadratic in 6a = a — (7^ ■ Indeed, calculating 6uj' = uj'[a] — u;'[crci] by 
the usual perturbation theory, we see these terms do not contribute in the first order 
because jo-^j = and / dz a'rjn = when nj^O. 

Minimizing the reduced effective energy ( p.6|) with respect to all possible configura- 
tions cr(x) one finds the mass of the soliton, M . A consistent approximation scheme for 
this procedure can be implemented as a power series in the coupling g that starts from 

M,i = minEei = EeiKi] = 0{g-^) . 

The next, Casimir, term modifies the mass by 

5M = E'cJa,,] = 0ig''). (3.16) 



Since, in general. 



SE'cJct] 



Sa 



g S^'kl^] 



2 \fl Sa 



^ , (3.17) 

"^cl/ regularized 



the Casimir term in E'^^ also changes the minimizing configuration, or the soliton form 
factor, by 0{g^). However, according to the previous paragraph, the non-local potential 
in eq. (|3.11|) does not affect the quantity (|3.16|) or the sum in eq. ( |3.17|) . We conclude 



that for the purpose of calculating the first quantum correction to the classical values 
of the soliton mass and form factor it is sufficient to sum up the oscillation frequencies 



3 



Since a;o = 0, it does not matter whether one includes it in the sum or leaves out. 



ojn about a spatially fixed configuration 0(x), omitting the non-local terms in eq. (|3.11|) , 
but throw away the lowest eigen-frequency, ujq . 

Nevertheless, eqs. ( |3.^ , p.llD allow one to compute the well defined one loop effective 
energy in the P = sector for the configurations not necessarily close to the solution of 
the classical equations of motion. This can be helpful in the search for non-perturbative 
radiatively generated solitons'^^^ , at least on the qualitative level. Therefore, in our 
calculational example in Section 6 we exploit the full form of eq. (|3.11| ) including all the 
non-local terms. 



4 Shape Function and the Soliton Form Factor 

In this short Section we demonstrate that, to our one-loop accuracy, the function cr(x) 
minimizing the functional -Egff ["^l ^^^ ^ physical interpretation of the soliton form-factor. 
Indeed, let us evaluate the matrix element of the operator 0(x) between the stable 
solitonic states with the mass M that differ by their total momenta: 

(P'|<^(x)|P) = e-(^'-^)-(P'|0(O)|P). (4.1) 

First, we present the state \P) as the soliton at rest, |O)0, boosted to the momentum P 
by the Lorentz boost operator Ap and recall the scalar nature of the field (/>(a;). We have 

(P'|0(O)|P) = (P'|0(O)Ap|O) = (P'|Ap0(O)|O) = (g|0(O)|O) = /(Q) , 

where 

|Q)=a1,|P')=A_p|P') 

describes a soliton carrying the momentum 
P'E — PF/ PP' 

Q = — j^^ = iP'-P)-^P'-P) + -. = P'-P + o{9^). 

For the very last estimate, we remember that M = 0{g^'^) and that we are interested in 
the momenta P and P' of the order of the inverse soliton size that is 0{g^) . 

Second, we substitute 0(x) by the decomposition ( \i.'S\) and Fourier transform the 
functions a{x) and ri{x): 

f{Q) ^ {Q\m\0)=9''{QWi-X)+gqkVki-Xm = 

= 9-'l^{Q\e''^'^HQ')+9qkUQ')]\0)- (4-2) 



The state |0) should not be confused with the true, perturbative, vacuum. 
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Third, let us notice that the state e ^^ ^\Q), being a momentum eigenstate because 

P (e-^«'^|Q)) = e-^«'^(P - Q')\Q) = (Q - Q') {e-"^'^\Q)) , 

is also almost an eigenstate of the Hamiltonian ( p.8|) : 

H {e-'^^'^lQ)) = e-'^^'^HlQ) + [H, e-*«'^]|Q) = M {e-'^'^\Q)) + 0(/)|g) . 

From these two observations we conclude that 

e-Q'^|g) = (l + 0(/)) \Q - Q') + 0{g')\{Q - Q'Y) , (4.3) 

where \{Q — Q')*) describes an excited soliton or a soliton-meson scattering state with 
the total momentum Q — Q' . The leading term in eq. ( ^73| ) is, of course, expected to be 
the same soliton boosted to the momentum Q — Q' for a heavy, non-relativistic particle 
described by the operators X and P. The higher order corrections should be anticipated 
as well because e*^"^ is not the complete boost operator in the field theory. 

Substituting the result (f4.3|) in eq. (f4.2|) and applying the normalization conditions^ 

(g-g'|o) = 27r5(g-g') , (4.4) 

(g-gigfe|o) = o, (4.5) 

one arrives at the formula 

f{Q)=g-'{^{Q) + 0{g'')] . (4.6) 



Summarizing eqs. ( |4.1| - |4.6| ), 



g{P'\4^{x)\P) = e-^(^'-^)-(a(P'-P) + 0(^3))= (4.7) 

'dye-^^'"-''^ya{x-y) + Oig') , (4.8) 



which is the ansatz (ii) of Goldstone and Jackiw^^l corrected for the one-loop quantum 
effects. In the perturbative expansion about the classical configuration, 

a{x) = aciix) + 5a{x) , (4.9) 



^ The normalization in eq. ( |4.4|) is taken to agree with the paper ||]. The second equation holds 
automatically for the forward matrix element, (0|(7fe|0), by our choice of the configuration a{x) that keeps 
the expectation values of g's at zero (see the text around eq. (^^)). For the off-forward matrix elements, 

{Q - Q'mo) = ^^ {Q - Q'\[p,qkm - , Q - ov . 
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we obtain the following linear equation for the leading order correction 5a that minimizes 
the sum Ed + -Ecas- 



"^^ +U'\ad{x)) 



dx 



gzri^^^))^^ (4,0) 



This is the formula for the one-loop correction to the soliton form-factor derived within 
the Kerman- Klein method'^! . The non-local terms in eq. (|3.11| ) take care of the (would 



be infinite) zero-mode term on the right hand side of eq. ( [4 .101 ) but give no other con- 



tribution in this order of the perturbative expansion ([4.9|) , as it was already stated at 
the end of Section 3. 



5 Regularization 

Now we return to the main line of our discussion and consider the reduced Casimir 



energy, currently written in the form (|3.8|) . The aim of this section is to convert this 
formal expression into an unambiguous and finite quantity calculable numerically with 
our computer. To avoid the many possible pitfalls in this step, we attempt to be very 
explicit in details. First of all, in order to regulate the oj'fjya] sum in eq. (p.8|), one 
should subtract the Casimir energy of the true, perturbative, vacuum. For the scalar 
field (t){x) living in a box of length L with fixed values at the boundaries, the normal 
oscillation modes ( "meson" excitations) about the topologically trivial minimum energy 
configuration (j){x) =g^^o'i (or o"2) have the frequency 

4°) = ^m^ + kl 



f^"^^^)+\-^^^] ' ^ = 0,1,.... (5.1) 



Each of the sums X) ^'k [^] ^"^^ Z) '^1°'* diverges badly and one should specify an ultra- 
violet cutoff prescription. It must be imposed consistently in both the trivial and the 
solitonic sectors^. This is not straightforward to achieve because of many incongruities 
between both the system description and its spectrum in the two sectors. 

Let us keep in mind that the ultra-violet divergences are unambiguously regularized 
by replacement of the continuous field 0(x) with a large but finite number of degrees 
of freedom having the same low-energy behavior, as provided, for example, by a lattice 
version of the theory. We imagine a continuous and physically realizable process that 
adiabatically transfers the regularized (quantum!) system from the vacuum state to 

^ It would be dangerous, and actually wrong in one space dimension, to cut off the oscillation modes 
at some fixed energy which is the same in both sectors because the number of modes below the cutoff 
in these sectors generally differs. 
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a(x) = tanh (x-a)/2 , a = -1.5, 0, 1.5, 5 




X -> 



Figure 1: Classical solutions of the Euler-Lagrange equation for a scalar field in the box 
< x < 10 and the (j)^ potential, eq. (^.j). The field satisfies the boundary conditions 
cr(0) = s£ [—1, 1] and cr(lO) ~ 1. The connection between the parameter s and the kink posi- 
tion, a, is given by eq. (|5.5|). 



a solitonic one by some external influence. As a specific example, we consider the 0^ 
modelQ 



U(a) 



a 



1 



(5.2) 



where we confine the scalar field 0(a;) to a finite box 0<a;<L,L^l and impose the 
boundary conditions 

4,{L)^g-^ 

that smoothly interpolate between the trivial and non-trivial topologies as the parameter 
s varies from 1 to —1. We would like to examine the behavior of the classical oscillation 
frequencies and the quantum Casimir energy in the ground state of this system as a func- 
tion of s. The classical time-independent field equation ( ^.15| ) with the potential ( |5.2| ) 
admits the solutions 



a[x] 



g<f){x) = tanh 



X 



Va = const 



(5.4) 



As the constant a in eq. ( p.4| ) varies from minus infinity to L/2 , cr{L) stays exponentially, 
i.e. up to 0(e~^), close to one, whereas a{0) goes from 1 to —1 + 0(e~^), as illustrated 
by Fig. 1. Therefore, up to unessential exponentially small corrections, we can satisfy 
the boundary conditions ( |5.3| ) by the solutions (§]^) where the parameter a is taken to 
be the following function of s: 

'1 -s' 



In 



1 + s 



(5.5) 



When the variation in s drives the system from the trivial sector to the topological 



one, the lowest of the normal frequencies ( |5.1| ) decreases from 
"^ The coefficients in eq. (^.2D are chosen so that ai^2 = Tl and U"{ai^2) 



LU, 



(0) 



l + (f) 



1 to 



1. 
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an infinitesimal value for L^oo. The zero point oscillations in the lowest normal mode 
become more and more broad until the quantum ground state of the system completely 
de- localizes. By this moment, the quasi-classical description of the corresponding normal 
coordinate breaks down, and it should be replaced by the collective quantization as 
s-^— 1. The other oscillation modes retain finite positive frequencies and their quasi- 
classical description holds all the time. Thus the first-order quantum correction to the 
soliton mass, M, may be determined as the sum of the shifts | Au;^ = ^ cj^ — | tu^ for 
the A; > 1 modes minus half of the frequency of the disappeared lowest vacuum mode. 
This is evident from the requirement of having an equal number of degrees of freedom, 
whether collectively quantized or not, in both sectors of the regularized theory. For 
notational convenience, let us define AtUg = Wq — cjg where co'g = corresponds to the 



zero eigenvalue of eq. ( |3.11| ), so that 



-| oo 

^1-loop = -Y.^^'k + C-t. . 
^ n=0 

In this equation "c.t." stands for the counter terms contribution due to one-loop renor- 
malization of parameters in the theory. We will return to it shortly. 

By now, regularization of E'Q^[a\ should become apparent. Indeed, from the phys- 
ical interpretation of the reduced effective energy, eq. (|3.5| ), we conclude that at its 
minimum, o"o, 

-| oo 

^Cas[^o] = Mi_ioop = - ^ Aa;^[ao] + c.t.[ao] . (5.6) 

^ n=0 

Since all the vacuum subtractions needed to regulate E'q^ [a] , together with the unphys- 
ical ci;Q[cr] =0, do not depend on the function (j{x) and we determined them for some 
particular configuration cro(a;), they simply carry over to an arbitrary ct{x) in the same 
topological sector as 

1 oo 

^CasM = :^EM.M+c.t.[a], (5.7) 

'^ n=0 

i^J,[a]^^M-^V ■ (5.8) 

Let us assume for simplicity that the Lagrangian ( p.l|) is invariant under the -^ —(f) 
transformation and take ai = —a2 for the a{x) boundary conditions (p.2|) . Then it is 
sufficient to consider only the configurations o"(x) that are odd functions of x so that the 



eigen-functions rin{x) in eq. ( |3.11| ) are either even or odd. Of course, now we imply that 
our regularization methods also respect parity, for example, the coordinate x varies from 
—L to L, where L :^ — . As x ^ — , rin{x) approaches the asymptotic form cos(A;„x + (5^) 
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or sin(A;„x + 5'_) for even or odd functions respectively with kn = J^jj'n^ — to? ■ The 
frequency shifts (^]8|) are related to the scattering phases 5^/c„) or 5'_{kn) as 




because kL = {k + Ak)L+5'^{k)^ provided k ^ 1/L. This part of the spectrum of eq. ( ^.11[ ) 
contributes to the sum ( p^ ) as 



j-i' (cont) 
-^Cas 



1 r I dkL\ ( duj 



2 J \ 71 ) \dk 
1 /■ kdk 



2tx J ,/W 



m 




(5.9) 



The spectrum of eq. (|3.11| ) will also have a few discrete modes with uj'^<m and the 



corresponding uj'^^ = m + 0{L ^). Summing up all the contributions, 



discrete 
spectrum 



The integral in eq. ( [5. 91 ) is a poor approximation when k < L^^ but because for such 
mode^ Auj' = 0{L~'^) and their number is 0{L^), we can safely set the lower limit of 
dk integration to zero, provided the integral converges in the infra-red. 

Since the scattering phases 6'^ fall off as 1/k for large momenta, the integral in 
eq. ( ^.10|) still diverges logarithmically on the upper limit. This ultra-violet divergence 



is removed by the usual parameter renormalization that manifests itself in the form of 
the counter term in eq. ( p.lO| ). Let us remember that counterterms are conventionally 
defined by a certain renormalization condition in the perturbative sector of the theory 
independently of the collective quantization procedure. By translational invariance, the 
corresponding term in the action, for instance. 



S^V oc J d'x [(f)\x) - 1\ (5.11) 

in the 0^ theory, can not depend on the collective coordinate X. The total momentum 
P is also ^f-suppressed in S'c.t.[0], as in eq. (|2.8| ), or does not appear at all in the absence 
of wave function renormalization that is the case in our example ( |5.11|) . Therefore, 



'S'c.t. = 'S'c.t. [?] simply modifies the potential V[q] in eq. (|2.8|), and for the purpose of 
obtaining the explicit form of the counterterm contribution to the reduced effective 
energy, we may ignore the reduction to the P = sector and consider the conventional 

^Unless the state becomes bound when it is counted separately. 
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* o 



p+k+1 

+ 



p+k 

Figure 2: The diagrammatic presentation of the one- loop contribution to the effective ac- 
tion, Fi^ioop- In the momentum representation, a hne carrying a momentum p stands for the 
usual free boson propagator -^^ — j, and the insertions are given by v{k) = J (Pxe^^^v{x) where 
k is the injected momentum. The integration goes over all the loop and the injected momenta, 
e.g. over p, k, and I for the third graph above. 



g V(x). 



effective energy E^s = T/T where T 
localized field configuration (j){x^ 

The one-loop part of F in the theory ( p.l|) is of the form 

'6^S 



is the usual effective action describing a spatially 



Ti- 



loop[ 



i In det 



; Tr In (-9- 



2 2 

— C — m 



Trln(-92-W" 

, v{g(j){x)) = U" {g(t){x)) - m^ 



(5.12) 



where again m^ stands for the value of U" at its minimum and describes the mass of 
the perturbative plane wave excitations over the true vacuum. Denoting the operator 
(—(9^ — m'^) by Aq, we have the following expansion in powers of the "interaction" f (0): 



Fi-ioopM = ^ Tr In Ao - « Tr (Aq^v) - % Tr (A^^'wA 



v\- ... 



(5.13) 



The first term in this expansion is a constant independent of the field configuration 
that is trivially removed by setting the vacuum energy to zero. The rest of the terms 
may be presented graphically by the sum of the diagrams shown on Fig. 2. The inclusion 
of Fc.tjf/'] = — 5'c.t.[</>] should remove all the divergences from the loop integrals on Fig. 2. 
For a scalar field in one dimension this can be achievedt^l by modifying the scattering 
phases in eq. (|5.1CI|) as 



b' 



6' -6. 



(1) 



Here 5+ and 5_ are the first Born approximations to the scattering phases in the 
continuous spectrum of the oscillation modes about a "nailed down" soliton proble: 



uj'^r]{x) 



5&^""W 



ri{x) 



-—^ + m^ + viaix)) 
dx'^ 



rj{x) 



(5.14) 



This trick works because the consecutive Born approximations for eq. (5.14) and the expan 



sion ( 5.13 ) are formulated as a power series of same quantity, v . Similar metliods can be carried out 
for sufficiently general theories of bosonic or fermionic fields in the realistic three space dimensions'^^! . 
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This operation is equivalent to adding a counter term that cancels the shift of the field 
vacuum expectation value in the perturbative sector due to one-loop corrections. The 
first Born approximations to the scattering phases in eq. ( p.l4| ) are easily found^'^^ as 



5+ (A;) = — -r dxv(x)cos kx 

k Jo 

S_ {k) = — — / dxv{x)sm kx 
k Jo 



and our final formula for the Casimir energy of the collectively quantized system reads : 

spectrum 

1 r-\-Qo 

where 4,(A;) = 6'_^{k) + 6'_{k) + - / dxv{x) . (5.16) 

rC Jo 

It should be reassuring that the non-local separable potential in eq. (|3.11|) does not affect 
the 1/k asymptotic behavior of S'_^{k) and the sum (|5.16|) falls off as 1/fc^ yielding a finite 



integral in eq. ([5.151). 



6 Numerical Computations 

As a demonstration of our method, we calculate the reduced effective energy for a set 
of field configurations in the one-dimensional 0^ example 

^ = ^(5.0)'-^M'-i)' + /:et . (6.1) 

We exploit two different families of trial configurations. The first one is the scaling 
transformation of the classical kink solution (|5T 



X 

ai(x) = g(t>{x) = tanh — . (6.2) 

2w 



The parameter w characterizes the width of the kink ( |6.2[ ) relative to its classical value. 
As it is seen from further calculations, the scaling transformation (|6^) possesses some 
unexpected peculiar properties. For this reason we also consider an alternative variation 
of the solution ( p.4| ) constructed to represent some "absolutely non-special" direction in 
the functional space of all admissible variations: 

X ( TIX \ 

o'liix) = tanh — h axexp . (6.3) 

The classical contribution to the effective energy, -Ecifo"], is trivially calculated from 
the formula ( W^ - Of course, the result should take the minimum, Mc\ = y2, at the 
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classical values of the variational parameters w = 1 or a = and grow quadratically with 
5uj and 5a, for example, 

E.[.,(x)] = ^ (w + -) = ^^r^-^ , 5w . w - 1 . (6.4) 

According to the previous sections, the reduced Casimir energy is determined by 
eqs. (|5.15| - |5.16| ) where the discrete oscillation mode frequencies, a;^, and the scattering 



phases in the continuous spectrum, 5'^{k), refer to the solutions of the auxiliary eq. (|3.11|) . 
The following Subsection 6.1 describes how to actually solve these non-local equations 
and extract the required quantities u'^ and 5'^{k). The reader who is not interested in 
these details should skip to Subsection 6.2 where we analyze our numerical results. 

6.1 Practical Implementation 

Substituting the explicit form of the potential W, eq. (|5.2| ), let us rewrite eq. ( p.ll|) as 



dy bi{y) r]{y) = [tv' - Ij i]{x) , (6.5) 



where we sum over the repeated index i = l,2 and define 

v{x)^U"{a{x))-l = ^-{a\x)-l) , 

ai{x) = ^=3Ux) , bi{y) = ^=a'{y), (6.6) 



a2ix) = ^=a'{x) , b2iy) = 

\/ J-rr 



/J-rr J-rr \/ J- rr 



with I, = J dzia'f, Ij = Jdzijy) . 

If the configuration a{x) is an odd function of x, which is true in our cases ( |6.2| , |6.3|) , 
the functions v{x), aj(x), and bi{x) are even and the modes {ri{x)} divide into the parity 
even (symmetric) and odd (antisymmetric) channels. The separable potential term in 
eq. ( |6.5D does not appear in the antisymmetric channel because the integral over dy 
trivially vanishes for the odd ri{y). Changing the lower integration limit to zero in the 
symmetric channel, we reformulate the problem (|6.5|) as: 

K^>V+ix) + 2ai{x) /o+°° dy bi{y)r]4y) = 

(symmetric channel) , (6.7) 
with ^7^(0) = 

Kuj'^-{x) =0 

(antisymmetric channel) (6.8) 
with 7]_ (0) = 



where K^^i is the local linear differential operator 

<P 
dx'^ 

First, let us determine the discrete spectrum, formed by the solutions with uj' <1 
that exponentially vanish as x ^ +00 . In the antisymmetric channel we simply integrate 
the differential equation ( |6.8|) from to +00 with initial conditions ^-(0) =0, ri'_{Q) = 1 
and search for u's for which ?7_ (x) falls off exponentially at large x . Of course, this has 
to be repeated for every configuration a{x). 

In the symmetric channel, eq. ( |6.7|) , we can also apply this shooting technique looking 
for the solutions in the form 

r]+{x) = ?7o(x) + air]i{x) , (6.9) 

where we sum over i = 1, 2 . In this ansatz rjo, rji, and 772 solve the ordinary differential 
equations 

K^'riq{x) = -ri'^{x)+v{x)r]q{x)-{uj''^~l)r]q{x) = \ ■ ,^ 

I aq{x) for q = t = 1, 2 

with the boundary conditions 

r/o(0) = l, r7;(0)=0; r],{0) , r][{0) = . 

Substituting the ansatz (|0|) into eq. ( |6.7|) we see that the latter holds provided the 
coefficients ai and 02 are determined from the following linear algebraic equation: 

tti + A + lij dj = , (6.10) 

where again i and j range over 1 and 2, and 

r+oo r+00 

A = 2 / dy bi{y) r]o{y) , 7^- = 2 dy hi{y) rjj{y) . 

Jo Jo 

The outcome of the described procedure applied to our first family of the trial configura- 
tions, eq. (|6.2|) , in both symmetric and antisymmetric channels is presented on Fig. 3 a) 
(solid lines). We would like to postpone its discussion until the following subsection. 

In the continuous spectrum of the symmetric channel we look for the solution of 
eq. (|6.7| ) in the form 

ry+(x) = e*^VW<^(a;) + e-'^'+^''^ ^* (x) 
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where ip{x) is a formal solution of eq. (|6.7] ) for x G [0, +00) such that 



(p{x) -^ e*^^ as X ^ +00 , 

and the scattering phase 5+(/c) is adjusted to satisfy the condition 

<(0) = . (6.11) 

In order to achieve better numerical precision, we trade the oscillating function (p{x) for 
a smoother one, c(x), as 



^, y^y<A^ J ^ 



(f{x) = c(x)e*''^ with c(+oo) = 1 and c{+oo) = 



(we remember that k = \Jlo''^ — 1 and c{x) depends on the parameter k as well). The 
function c{x) also satisfies a linear differential equation with a separable potential. 



- d'{x) - 1%kd{x) + v{x)c{x) + 2e-*'=^ai(a;) / dy bi{y)e"'yc{y) = , 

Jo 

and is determined, just as in the case of eq. (|6.7[), with the ansatz 

c(x) = Co(x) + ajCj(a;) , (6-12) 

where 



c^{x) -2ikCq{x)+v{x)cq{x) , _ ,__, ^.^^^ 



for g = , 

for g = 1, 2 ; 



co(+oo) = 1 , Co(+oo) = ; q(+oo) , c-(+oo) = . 

Again, we have the algebraic equation a, + /3j + ^ij aj = for the coefficients ai and a2 
where 

A = 2 / dy h{yy^y c,{y) , %• = 2 / dy 6,(i/)e^'=^ c, (y) . 



Having calculated the function c(x), we find the phase 6'^{k) from the condition ( |6.11| ) 



as 



5; = -arg[A;c(0)-ic'(0)] , 

where the branch of the argument functions is unambiguously specified by the require- 
ment limfe^oo 5^ = . 
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Figure 3: The normal oscillation frequencies in the discrete spectrum, (a), and the reduced 
Casimir energy together with the contribution to it from the continuous spectrum, (b), versus 
the soliton width, w, in the ansatz aj = tanh ^ . The solid lines on the plots describe cv'^, 
-^Cas ' ^^^ -^Cas ^^ ^^^ theory reduced to the P = sector. The dashed lines are the results 
in a conventional analysis of small oscillations about a kink with a defi- 



for uJr, and E, 



(cont) 
Cas 



nite position in space. The only drastic difference between the two approaches occur in the 
translational mode, for which lo'q = 0, but ujq oc \/1 — w and is imaginary when w > 1 . 



The antisymmetric channel is trivial: The continuous solutions of eq. ( |6.8|) are given 



by 



where 



and 



rj_{x) = e*^'-(*^V(a;) - e-'^'-^''^^*{x) ( such that 77_(0) = ) 



(p{x) = c{x)e 



ikx 



with c(+oo) = 1 , c'(+oo) = 



— c"{x) — 2ikc{x) +f (x)c(x) = . 

From the requirement r]^ (0) = we obtain 

6'_ = — arg c(0) 

in the argument branch such that linifc^oo S'_ =0 . Calculating the scattering phases in 
both channels and regularizing them by removing the first Born approximation, eq. ( |5.16| ) 
we obtain the continuous spectrum contribution to the Casimir energy. 



6.2 Discussion of the Results 

The spectrum of small oscillations about the classical kink solution o"ci = tanh | consists 
of the continuous spectrum with u = \/l + k'^, a parity odd discrete state at a;_ = ^, 
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and the even ti;o = translational mode. Besides, the parity even amphtude 77+(x) cor- 
responding to ujj^ = 1 approaches a constant value as a; ^ oo. Infinitesimally lowering 
the frequency with an arbitrarily small variation of the soliton shape, one would send 
this state into the discrete spectrum. It is often convenient to refer such a state as a 
half-bound one. 

A small variation of the shape function induces a shift in the oscillation frequencies 
and the total Casimir energy. In general, the shifts are expected to be linear in the 
variational parameters, (5w = w — 1 or a in our examples ( |6.2| ) or ( |ti.3D , of course, up to 
less important in the perturbative regime higher power corrections. The non-vanishing of 
-^^ is clearly visible for the odd discrete state on Fig. 3 a) and for the continuous 

contribution to the Casimir energy on Fig. 3 b), where the ansatz ( |6.2| ) was considered. 
The dashed lines on this Figure present the same quantities computed for a localized kink 
configuration'^'^*'', when the non-local terms in eq. (|3.11| ) are absent. The discrepancy 
between the two methods for the non-zero modes is barely observable, demonstrating 
the weakness of the collective quantization influence on these modes in the perturbative 
regime. 

As expected, there always exists a normalizable zero-frequency solution in the sym- 
metric channel of the reduced problem (|3.11| ). For a localized configuration, the "zero"- 



mode frequency would vary as ujq oc yjl — w and become imaginary as w exceeds 1 . This 
square root singularity and the imaginary contribution to the conventional effective en- 
ergy has nothing to do with the soliton instability (it is stable!). It is rather due to 
inadequacy of the standard effective energy formalism in describing the true solitonic 
ground state that has no specified position in space. 

The wide solid line on Fig. 3 b) is our final result for E'q.-^ verses the kink width. The 
first perturbative correction to the soliton mass'^'^^'^"', 

AAfi_,„,, = 41(0) = -1= - A ^ -0.333127 , (6.13) 

is numerically found to be —0.333124, which, given our numerical accuracy, is in com- 
plete agreement with the analytical result (|6.13|) . Surprisingly, the plot on Fig. 3 b) 



shows no apparent dependence of -E^as ^^ ^^^ width w. Indeed, in Table 1 we give the 
deviations of numerically computed contributions to -Ecas froni their classical values for 
the antisymmetric discrete mode uj'_{w), which coincides with ^^-(w), and for the sum 
of the continuous modes, -Ecas" i ^^e latter being separated into -EcaT obtained for a 
localized soliton and the correction to it due to collective quantization. As one can see 
from the table, the linear variation in the discrete spectrum exactly cancels by the local 
part of the continuous spectrum and one is left with the quadratic in 5w dependence, 
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Table 1: The change in the reduced Casimir energy (the last line) against the soli- 
ton width variation and the main three contributions to it. The results demonstrate 
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j{cont) j-^{cont) 



and 



{Eq^ — Eq°^ ) oc (w — 1)^ . The last two values for -E^as' marked by an asterisk, should 
also include the symmetric discrete mode with the frequency very close to one. 
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Table 2: The numerical results for the second family of trial configurations, {cr//(a;)}. 



which is largely produced by the non-local terms in eq. ( |3.11D . Thus our calculations find 
no one-loop quantum correction to the soliton width, w. This result comes unexpected, 
remembering that cr(x) should receive non-vanishing one-loop correction according to 
eq. ( [4.10 ). We do find such a correction when we consider another family of trial con- 



figurations, {o"7/(x)} in eq. (O). The numerical results for this choice are summarized 
in Table 2. We obtain that 



dE'r 



Cas 



da 



dE, 



Cas 



«=0 



da 



-1.5 ±0.1) X 10" 



Q = 



resulting in 



6ai- 



loop 



(^) 






/• (2.2 ±0.2) X 10 



-2 



a=0 



We see that the first quantum correction to the shape parameter a is small until g^ <^ 10^ . 
However, the contribution of the non-local terms becomes significant at a much smaller 
couphng. From our calculations. 
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,/{cont) 
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4.2 ±0.1 



a 



(6.14) 



whereas g^i^ ^0.6956/^^ . The effect of the quantum correction ( 6.14 ) becomes com- 



parable with the classical energy at g'^ 



4.2 



■0.17. 
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7 Conclusion 

To summarize, we adapted the conventional effective energy variational techniques to 
the situation where the Hamiltonian of the system possesses cyclic variables and the 
standard semi-classical methods should incorporate their collective quantization. In 
our formalism one may retain the description of the system in terms of a continuous 
field, (l){x), by introducing an external source with a special non-linear non-local cou- 
pling. The source term in eq. ( 3.14|) insures the coupling to all the normal degrees of 



freedom excluding the cyclic variables. The effective energy functional reduced to the 
P = sector of solitonic states, E'^ff[a{x)], is given to first order in h by the classical en- 
ergy of the configuration 0(x) =g~^(j{x) and the properly regularized reduced Casimir 
energy, eq. ( ^.15| ). The mass of the soliton in the ground state equals the minimum 
of £'gfj[cr], and the configuration a{x) at which this minimum is achieved describes the 
soliton form- factor, eq. ([4.8|). 

The effect of collective quantization technically amounts to the introduction of the 
non-local separable potential]^ in the eigen-mode equation ( |3.11|) . However, in order 



to obtain the first quantum corrections to the classical values of the soliton mass and 
form-factor, one could ignore these non-local terms but follow the rule that the lowest 
oscillation eigen-frequency must be excluded from the Casimir energy. The correction to 
the reduced energy due to collective quantization grows quadratically with the deviation 
from the classical field configuration. To get a possible insight into the physics in the 
non-perturbative regime, one may still consider £'gg[o"(x)] computed only to the one-loop 
order. This can be done numerically for an arbitrary field configuration as described in 
Section 5. 

It would be interesting to apply this method to investigation of the role of collective 
quantization in deeply non-perturbative regime. This could be necessary, for example, 
for a consistent treatment of the models^^^' in which a non-classical field configuration is 
stabilized by an order h contribution to the soliton energy from fermionic fields coupled 
to the scalar "Higgs" . An extension of the method that would incorporate the fermionic 
fields as well, presents another challenge by itself. 



^° Some aspects of the Casimir energy for systems described by separable potentials has been con- 
sidered by Jaffe and Williamsonl^^l . 
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